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We consider the existence of bound systems consisting of two quarks and two 
' antiquarks (q 2 q 2 ) within the framework of a constituent quark model. The 

underlying quark dynamics is described by a linear confinement potential and 
an effective q 2 q 2 interaction which has its origin in instanton effects of QCD. 
in . We calculate the spectra and examine the internal structure of the states 

^ . found. 

in 
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in ■ I. INTRODUCTION 
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Abstract 



Most hadrons found experimentally have been identified as mesons (qq) or baryons (q 3 ) 
and their masses and decays can be described fairly well in various quark models. However, 
there are still a few so-called exotic mesons which can not be interpreted as qq states. For the 
lightest exotic mesons /o(975) and ao(980) an underlying q 2 q 2 structure (or KK molecule), 
besides other interpretations (e.g. has been discussed but could not been proved 

yet unambiguously. Recent experiments and data has renewed the interest in these exotic 
states 

We will discuss the existence of bound q 2 q 2 states within a non-relativistic constituent 



quark model, which was successfully used to describe meson and baryon mass spectra flO 
and decays. We can compare the states of this four particle problem with uncorrelated two- 
meson states and thus obtain a first clue whether the calculated states in reality would be 
bound or dissociated. In addition, the model is just simple enough to make the calculation 
of this real four particle problem possible. 

We consider the three lightest quark flavours u, d and s, where the constituent quark 
masses of n and s (n stands for a u or d) are assumed to differ. The potential part of the 



*Submitted to Z. Phys. A 

^E-Mail: beinker@ptprslO.phy.tu-dresden.de 



1 



Hamiltonian consists of a sum of two-body forces. This contains a long range, colourdepen- 
dent confinement part and a short range instanton induced interaction, which is generalized 
to three flavours [[K] . 



The hamiltoninan was calculated in a full basis of solutions of the 3 x 3-dimensional 
harmonic oscillator. Solving the eigenvalue problem and by a variational principle yield 
the mass spectra and the corresponding eigenstates for further examination of the internal 
structure of the states. 



II. THE MODEL 

For calculating the q 2 q 2 states, we used a Hamiltonian of the form 

H = M + K + V con{ + W (1) 

where 

M = Y,m i 

i 

is the sum of the four constituent quark masses and 

p 2 p2 



2m,- 2M 



the kinetic energy for the relative motion. 

As confinement potential Vconf, we choosed a linear rising potential with the relative 
quark- (anti) quark distance: 

i<j 



The offset and the slope b^ are different for qq and qq confinement. As noted in JTT 
a qq and b qq are not independent. The dependence of b qq from b qq is given by geometrical 



reasons [13]. The following relations should hold: 



b qq = 0.5493-\ g - 

Although the parameters were independently varied, the above relations in the present 
approach are fullfilled with great accuracy. They are given in table ||[ We consider two 
different kind of colour dependencies Fc of V^ nf, which in the meson or baryon case are 
identical. In the first case, Fc is identical to the colour singlett projector P c [[|||: 

Konf ( r ij) = Pc( a ij + ^ij r ij) 
i<j 

whereas in the other case we have: 
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i<j 

with Aj is a Gell-Man matrix. The V^, nf potential acts only on colour singletts. In the case 
of q 2 q 2 only the four qq pairs will be affected. Since it is an purely attractive potential 
two distinct qq colour singlett states (which we call quasi-mesons) will always have a net 
attraction, because also qq pairs with q and q from different quasi-mesons have always colour 
singlett contributions yielding unwanted long range forces. 

The V^ nf potential has a better behaviour for R — > oo, R is the distance of the two 
quasi-mesons. As was shown [13 , it goes with R~ 3 as R grows large which is a consequence 
of the repulsive part of V^ nf , which is an effect without experimental evidence. In addition 
V^ nf allows a better comparison with earlier results of other authors @J|. 

Finally, 

W = Hq 1 q 2 + Hg s g 4 + H qi q 3 + Hq 2 q A + H qi q 4 + Hq 2 q A 

is the residual quark- (anti) quark interaction. The q\q^ term is of the form 

H qiq2 = -~9(P S ^Pe C + 2P**P? 

where the flavour matrix g is shown in table |H], P 5 are spin projectors on the relative qq 
spin state and P c denotes a colour projector, r is the relative distance of the two quarks. 
The qq term is of the same form, if one changes all colour to anticolours and all flavour to 
their antiflavour indices. All qq terms are of the form 

H qq = g (jp^Pf + P S =° QP 8 C + 8Pf )) 5 3 (r) 

where the flavour depending g is given in table |TT| (for the calculation of this residual 
interaction from an effective Lagrangian see |10|| ). 

The matrices g and g include the total flavour dependence of W. As one can see from 
table U, 9 contains a projector on flavour antisymmetric states, while g generates a flavour 
mixing for T = states, which allows to describe the rj-rf splitting in the meson case, and 
lowers the energy of the isovector state. 

The more complicated structure of W for q 2 q 2 in comparison to the meson or baryon 
case is a consequence of the fact that there exist two distinct colour singletts for q 2 q 2 and 
consequently the terms including Pg or P$ do not vanish. The residual interaction is a pure 
contact force and as such leads to an unbound Hamiltonian. Therefore, we have chosen to 
regularize the interaction by replacing the 5 distribution by a Gaussian 

A" 7T2 

where A can be interpreted as an effective range of the pairing force. We note that a 
derivation of the 't Hooft force going beyond the one loop approximation which was used, 
is expected to such a finite range |TU| . 



For the calculation of the q 2 q 2 mass spectra and the eigenstates, we proceed as follows: 
The matrix elements of the Hamiltonian ([!]) are calculated within a spin-flavour SU (6), 0(3) 
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oscillator basis for each relative coordinate comprising N = 8 oscillator excitations in total. 
We should note that this size of the model space is still too small for reaching convergence 
of the spectra. 

Calculations with N > 8 oscillator excitation were limited by the calculational effort 
However, the comparison with meson calculations shows that the lowest eigenvalues are 
calculated with an accuracy of a few MeV. 

We obtained the q 2 q 2 spectra by application of the variational principle, i.e. diagonal- 
isation of the Hamiltonian in the oscillator basis and minimization of the lowest energy 
eigenvalue with respect to the oscillator length parameter. This is consistent to the earlier 
meson and hadron calculations in |]I0 |. The parameters entering the present calculation are 
given in table |I[ These parameter were fitted to the meson and baryon spectra. 



III. THE SPECTRA 

Fig. [I] shows the mass spectra in the most interesting scalar isoscalar and isovector 
channel. The columns for each channel are devided into four subcolumns, which show from 
left to right the sum of the experimental masses of two mesons, the sum of two meson masses 
calculated with the the same model in the meson sector, the calculated q 2 q 2 spectra with 
V^ nf and, finally, the spectra calculated with V^ nf . In addition, the very left subcolumns 
shows the masses of /o(975) and ao(980), the most popular q 2 q 2 candidates. 

The q 2 q 2 eigenvalues are denoted by n 2S+1 L pi2PS4 where S is the total spin, L the total 
angular momentum and P12 resp. ^34 the exchange symmetric of q\ with q2 and 53 with q^ 
respectively, n counts all states with identical S, L, p 12 and p 34 . 

The two meson-state masses, which are compared with the q 2 q 2 spectra, are choosen to 
have the same observable quantum numbers as the q 2 q 2 states with the two mesons in a 
relative s-wave. A main diffiulty of our q 2 q 2 model is, that we can not describe two free 
mesons in a finite model space with a finite number of oscillator excitations. Otherwise, 
it would be impossible to calculate the Hamiltonian. It is therefore neccessary to examine 
the internal structure of the q 2 q 2 states in detail to differentiate real bound q 2 q 2 states from 
quasi free two meson states. 

As one can see from fig. |l|, the q 2 q 2 spectra describe in good approximation the lowest two 
meson mass sum for the spectra calculated with the V^, nf potential. The lowest state O ^ 
lies under the m n +m n treshold which is due to the strong long range meson-meson attraction 
from the V^, nf potential. Since we consider this long range attraction to be unphysical, this 
does not imply that O ^ is in reality a bound q 2 q 2 state. The eigenstates calculated with 
the V^, nf , especially the lowest one, have higher energies due to the repulsive part of this 
potential. 

As a consequence of the discretization of the two meson continuum in our variational 
approach, both spectra contain more eigenstates in the energy region shown than the sum of 
the masses of two mesons (with the two mesons in a s-wave) would predict. In other words, 
the spectra alone do not allow to draw any conclusions wether there exist any "real" q 2 q 2 
bound states within this model or not. A closer examination of the internal structure of the 
calculated eigenstates is therefore neccessary. 
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IV. INTERNAL STRUCTURE 



To obtain such a closer view of the interpretation of the internal structure of the cal- 
culated q 2 q 2 states, we calculated various probability densities. Let us consider the lowest 
eigenstate 0°>S . (Unless stated otherwise we refer to the calculation with the V^ on{ po- 
tential.) Fig. shows the densities for the various relative distances with N max = 2 and 8. 
These densities were calculated by summing over all internal quantum numbers and inte- 
grating the q 2 q 2 probability density ^^(ui, u 2 , u 3 )^ a (ui, u 2 , u 3 ) over all relative coordinates 
except the relative distance under investigation, with U; L denoting a relative coordinate (e.g. 
the g-g-distance) and a a multiindex for all internal quantum numbers. 

The qq distribution is calculated as: 

Pqq{ u ) = / d U q2 ^q i d U qi q 3y q 2 q i \ty a (u, U q2 ^ 4 , U qi q 3> q 2 q 4 ) | , (2) 

a J 

since all four pairs have the same distribution. It is observed that all densities in fig. |2| 
are spreading with growing N max . While the q-q distance (or inner quasi-meson radius) 
distribution keeps its form, the qq-qq (or quasi-meson) distance distribution is mostly spread 
and its form seems to go over in a uniform distribution with rising N max , only slightly 
modified by the long range part of the V^, nf potential, as mentioned above. 

For higher excited states the qq-qq distribution gets more complex and broader while the 
q-q distribution again ist mostly localised, and the same conclusions as above can be drawn. 
Comparison with the distribution of 7r-meson radius within the same model shows that it 
has nearly the same form as the q-q distribution of the lowest q 2 q 2 state 0°S (see fig. §). 

At last, we considered the distribution p(z) of a single quark on the z-axis given by: 

p(z) =J2 J d 3 xid 3 x 2 d 3 x 2 d 3 x 2 \^ a (xi, x 2 , x 3 , x A )\ 2 5 rrnx^j 8(xi - x) (3) 

where £j are the absolute coordinates of the ( ant i-) quarks. The first ^-function ensures that 
the center of mass is identical with the center of the coordinate system, and 5(xi — x) cancels 
the integration over x\\ x is defined by: 



x = . (4) 



Surprisingly, the p(z) depends only very weakly on N max , so we can neglect this dependence 
for a qualitive analysis. The distribution has a quite sharp maximum at z = and two 
significant shoulders at z — ±lfm. If we make the assumption, that we have a superposition 
of two Gaussian-like quasi-meson distributions with their maxima seperated by a distance 
2a, we would expect a distribution of the form 

e -(x+af + e -{x-a? = e -^-a* (^+2ax + e ~2a^ = ^ cosh ( 2a x). (5) 

So we would expect a Gaussian which is modulated by a cosh-function. Now it is possible 
to write p(z) in the form 
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z 

p(z) = Y^ a i( z ) e ^ 

i 



(6) 



where the width hi is a flavour dependent constant. The coefficient functions di(z) corre- 
sponding to the Gaussian expansion @ are shown in the fig. |3]. 

It is clear to see, that for z — > ±00 the coefficients which give the important contributions 
to the sum behave much like cosh functions. This gives evidence that the 0°S state in fact 
describes two free quasi-mesons and is not a bound q 2 q 2 system. Although the contributions 
of the different flavour components change, this picture is qualitatively still correct if one 
considers the next higher excited states, also for the V^ nf . That means, all states could be 
most-likely described as free quasi-mesons. 

The little local maximum of the n 2 n 2 coefficient in fig. [3] may indicate that there exist a 
meson-meson attraction for small distances, which is, however, too weak to lead to a bound 
system. 

V. CONCLUSION 

We have extended a model, which is well established in the meson and baryon sector, to 
q 2 q 2 systems, without any change in the interaction. 

The interpretation of the results was not easy, because of the difficulty to distinguish 
q 2 q 2 from bound states and quasi-free mesons. On the basis of the internal structure, i.e. the 
density for a single quark and its representation as a sum of Gaussians, we believe, however, 
to have found a good method to identify possible q 2 q 2 states in spite of the background of 
quasi-free meson states. 

On the basis of results presented in this contribution, the only plausible interpretation is 
that all states we calculated correspond to quasi-free mesons. Especially, we found no hint 
that there exist very tightly bound states of two quarks and two antiquarks. However, this 
does not yet rule out, that the Jo (975) and ao(980) mesons are KK molecular states, bound 
by long range forces, as provided e.g. by meson exchange 0. 
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FIG. 1. Spectra for scalar-isoscalar and vector-scalar q 2 q 2 -systems. The subcolumns show 
from left to right the sum of the experimental masses of two mesons, the sum of two meson masses 
calculated with the the same model in the single meson sector, the calculated q 2 q 2 spectra with 
V^ on{ and, finally, the spectra calculated with V^ nf 
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FIG. 2. Density distributions in the channel 0(S + ) in dependence of N max . The distributions 
belong to the ground state. A Q denotes here an anti-quark q. For comparison the q-q distribution 
of the 7r-meson calculated in the same framework is shown. 
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FIG. 3. Coefficient functions belonging to the Gaussian expansion of p(z) in (^) for the 
different flavour components with N max = 2. n denotes an u or d quark, a capital letter the 
corresponding anti-quark. 
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TABLES 



TABLE I. The parameters of the Hamiltonian, fitted to the mass spectra of the light mesons 



and baryons [10 



m n 


300 


Mev 


m s 


540 


Mev 




-892 


Mev 


bqq 


850 


Mevfm" 1 


0>qq 


-511 


Mev 


b qq 


467 Mevfnr 1 


9 


122 


Mev fm a 


9' 


82 


Mev fm 3 


A 


0.37 


fm 



TABLE II. The flavour matrix g with g' := |ff c ff(n) and g := fg c ff( 5 )- 





ud du 


us 


ds 


sd 


su 


uu dd ss 


ud 
du 


9 ~9 
-9 9 








us 




9' 








-9' 




ds 








9' 


-sf 








sd 





-9' 


9' 





su 




-9' 








sf 




uu 














dd 

















ss 















TABLE III. The flavour dependent matrix g with g and g' the same as in table p. 





ud du 


us 


ds 


sd 


su 


uu 


dd 


ss 


ud 


-9 


















du 
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